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I review recent lattice calculations of hyperon physics, including hyperon spectroscopy, axial coupling constants, 
form factors and semileptonic decays. 



1. Introduction 

The hyperons are extremely interesting because 
they provide an ideal system in which to study 
SU(3) flavor symmetry breaking by replacement 
of up or down quarks in nucleons by strange ones. 
Hyperon semileptonic decays provide an addi- 
tional way of extracting the CKM matrix element 
V us and offer unique opportunities to understand 
baryon structure and decay mechanisms. How- 
ever, since hyperons decay in less than a nanosec- 
ond under weak interactions, their experimental 
study is not as easy as for the nucleons, and thus 
hyperon properties are not as well determined. 
Recent successes of lattice QCD in computing nu- 
cleon structure |H2l3j provide assurance that lat- 
tice QCD can reliably predict the properties of 
hyperons as well. Knowledge of these properties 
can be very valuable in understanding hypernu- 
clear physics, the physics of neutron stars and the 
structure of nucleons. 

In this proceeding, I review the latest progress 
on the hyperon spectroscopy, axial coupling con- 
stants, electric charge radii, magnetic moments 
and semileptonic decay form factors using a 
mixed action, in which the sea and valence 
fermions use different lattice actions. Lattice 
calculations of hyperon scattering lengths are 
not included in this proceeding; for more de- 
tails, please see Ref. [I]. In our case, the sea 
fermions are 2+1 flavors of staggered fermions (in 
configuration ensembles generated by the MILC 
collaboration^ ), and the valence fermions are 
domain- wall fermions (DWF). The pion mass 
ranges from 300 to 700 MeV in a lattice box of size 
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2.6 fm. The gauge fields are hypercubic-smeared 
and the source field is Gaussian-smeared to im- 
prove the signal. Details on the configurations 
can be found in Ref. [6]. For three-point func- 
tions, the source-sink separation is fixed at 10 
time units. 

2. Spectroscopy 

On the lattice, continuum SO (3) rotational 
symmetry is broken to the more restricted sym- 
metry of the cubic group (also known as the oc- 
tahedral group). Hadronic states at rest are clas- 
sified according to irreducible representations (ir- 
reps) of the cubic group. Since cubic symmetry 
is respected by the lattice action, an operator be- 
longing to a particular irrep will not mix with 
states in other irreps. The most general baryon 
spectrum from a lattice calculation is given by 
a small number of irreps of the double-cover of 
the cubic group: G lg/u , H g/u and G 2g / U , where g 
(German: gerade) and u (ungerade) denote posi- 
tive and negative parity, respectively. 

We follow the technique introduced in Ref. [7] 
to construct all the possible baryon interpolating 
operators that can be formed from local or quasi- 
local u/d and s quark fields; the G^g/u irreps are 
excluded, since they require non-(quasi-)local op- 
erators. In this paper, we present the calcula- 
tions of the masses of the lowest-lying states in 
the Gi g / U and H g / U representations. 

The non-locality in four dimensions of our va- 
lence DWF action manifests in oscillations of the 
two-point effective mass close to the source. As 
a result, a phenomenological form for fitting such 
data was proposed in Ref. [5] , employing both an 
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oscillating contribution describing the non-local 
lattice artifacts and two positive-definite contri- 
butions: 

C(t) = Aoe-W-*™) + A l e- M ^ t ~ t ^ 
+ A^i-lfe-^-^. (1) 

The first excited-state mass Mi is included to ex- 
tract a better ground-state mass M ; M osc is a 
non-physical oscillating term due to lattice ar- 
tifacts. We use two-point correlators with the 
same smearing and interpolating operators at 
both source and sink and select fitted results with 
varying fit ranges, optimized for quality of fit. A 
standard jackknife analysis is employed here. 

Here we will focus on orbitally excited hyperon 
resonances. The ground states of the octet and 
decuplet and the SU(3) Gell-Mann-Okubo mass 
relation can be found in Ref. 6 . There is in 
this reference an extensive description of baryon- 
mass extrapolation (to the physical pion mass) 
using continuum and mixed-action heavy-baryon 
chiral perturbation theory for two and three dy- 
namical flavors. However, no chiral perturbation 
theory results have been published for orbitally 
excited hyperon resonances; therefore, we apply 
naive linear extrapolations in terms of M%. Fig- 
ure [1] summarizes our results for the lowest- lying 
A and f2 Gi g / U (upward/downward-pointing tri- 
angles) and H g / U (diamonds and squares) (since 
these have less overlap with the data in in Ref. [B] ) 
and their chiral extrapolations. The leftmost 
points are extrapolated masses at the physical 
pion mass, and the horizontal bars are the ex- 
perimental masses (if they are known). 

We summarize the lattice hyperon-mass calcu- 
lations for the E, A, 3 and fl in Figure[2j divided 
into vertical columns according to their discrete 
lattice spin-parity irreps (Gi g / U and H g / U ), along 
with experimental results by subduction of con- 
tinuum J p quantum numbers onto lattice irreps. 
G\ ground states only overlap with spin- 1/2. The 
spin identification for H can be a bit trickier, 
since this irrep could match either to spin-3/2 or 
5/2 ground states. We simply select the lowest- 
lying of 3/2 or 5/2 indicated in the PDG, so it 
could be either depending on which one is the 
ground state for a particular baryon flavor. (We 
note below if the lowest H is not spin-3/2.) 
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Figure 1. The squared pion- mass dependence of 
A and ^-flavor baryons with their extrapolated 
values. The bar at the left indicates the experi- 
mental values for the corresponding spin. 



Although our naive extrapolation neglects con- 
tributions at next-to-leading-order chiral pertur- 
bation theory, several interesting patterns are 
seen. The better-known A (with H g being spin- 
5/2) and E spectra match up with our calcula- 
tions well. Giu 2 lines up well with the S(1690), 
indicating that the spin-parity of this resonance 
could be 1/2 - ; this agrees with SLAC's recent 
spin measurement. The spin assignments for the 
£1 channel are the least known; from our mass 
pattern, we predict that Sl(2250) is likely to be 
3/2~ (although we cannot rule out the possibil- 
ity of 5/2-), and 0(2380) and 0(2470) are likely 
spin 1/2 - and l/2 + respectively. 

To successfully extract any reliable radially ex- 
cited hyperon states, we would need finer lattice 
spacing, which would require massive computa- 
tional resources to achieve. In Euclidean space 
the excited signals exponentially decay faster 
than the ground state. One possible solution is to 
use an anisotropic lattice, where the temporal lat- 
tice spacing is made finer than the spatial ones to 
reduce overall costs. A full-QCD calculation us- 
ing anisotropic lattices is underway, and we would 
expect new results on these excited states within 
the next couple of years. 
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as 350 MeV. Figure [4] shows our lattice data 
as a function of f-w) 2 with the correspond- 
ing chiral extrapolation; the band shows the 
jackknife uncertainty. We conclude that qa — 
1.18(4) stat (6) syst , gss = 0.450(21) stat (27) syst 



and 



-0.277(15) stat (19) 



syst- 



In addi- 



tion, the SU(3) axial coupling constants are es- 
timated to be D = 0.715(6) st at(29) syst and F = 
0.453(5) st at(19)syst- The axial charge couplings of 
£ and 2 baryons are predicted with significantly 
smaller errors than estimated in the past. 

We also study the SU(3) symmetry breaking in 
the axial couplings through the quantity ^su(3): 



^SU(3) 



g A -2.0xg ss +g E 



(2) 



Figure 2. A summary of our hyperon spectrum 
compared to experiment. Our data are the short 
bars with large labels on the right. The experi- 
mental states are long bars with small labels on 
the left. 



3. Axial Coupling Constants 

The hyperon axial couplings are important pa- 
rameters entering the low-energy effective field 
theory description of the octet baryons. At the 
leading order of SU(3) heavy-baryon chiral per- 
turbation theory, these coupling constants are lin- 
ear combinations of the universal coupling con- 
stants D and F, which enter the chiral expansion 
of every baryonic quantity including masses and 
scattering lengths. These coupling constants are 
needed in the effective field theory description of 
both the non-leptonic decays of hyperons, and the 
hyperon-nucleon and hyperon-hyperon scattering 
phase shifts [9]. Hyperon-nucleon and hyperon- 
hyperon interactions are essential in understand- 
ing the physics of neutron stars, where hyperon 
and kaon production may soften the equation of 
state of dense hadronic matter. 

We have calculated the axial coupling con- 
stants for X and 3 strange baryons using lat- 
tice QCD for the first time. We have done 
the calculation using 2+1-flavor staggered dy- 
namical configurations with pion mass as light 



where x is (Af| - M%)/(4nf%). Figure [3] shows 
^su(3) as a function of x. Note that the value 
increases monotonically as we go to lighter pion 
masses. Our lattice data suggest that a <5su(3) ~ 
x 2 dependence is strongly preferred, as the plot of 
6su(3)/ x2 versus x in Figure [3] also demonstrates. 
A quadratic extrapolation to the physical point 
gives 0.227(38), telling us that SU(3) breaking is 
roughly 20% at the physical point, where x = 
0.332 using the PDG values^ for M„+, M K+ 
and fn+ ■ We compare the result of heavy-baryon 
SU(3) chiral perturbation theory [TT] for <5su(3) as 
a function of x, and we find that the coefficient 
of the linear term in Eq. [2] does not vanish. This 
implies that an accidental cancellation of the low- 
energy constants is responsible for this behavior. 

4. Form Factors 

The study of the hadron electromagnetic form 
factors reveals information important to our un- 
derstanding of hadronic structure. The electro- 
magnetic form factors of an octet baryon B can 
be written as 



(B\VJB)=u B 



Jtj.F^q 2 ) + a^qy 



F 2 (q 2 ) 



2M B 



u B 



from Lorentz symmetry and vector-current con- 
servation. Fi and F 2 are the Dirac and Pauli 
form factors. Another common form-factor defi- 
nition, widely used in experiments, are the Sachs 
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Figure 3. (Top) The SU(3) symmetry breaking 
measure 5su(3) • The circles are the measured val- 
ues at each pion mass, the square is the extrapo- 
lated value at the physical point, and the shaded 
region is the quadratic extrapolation and its er- 
ror band. (Bottom) Sgu^/x 2 plot. Symbols as 
above, but the band is a constant fit. 



form factors; these can be related to the Dirac 
and Pauli form factors through 



T 



G E (q 2 ) = F.iq 2 ) 
G M (q 2 ) = F^q 2 ) + F 2 (q 2 ) 



AM 2 



(3) 
(4) 



In this work, we will concentrate on Sachs form 
factors. Note that here we only calculate the 
"connected" diagram, which means the inserted 
quark current is contracted with the valence 
quarks in the baryon interpolating fields. 

On the lattice, we calculate the quark- 
component inserted current, V M = <n M <z, with 
q = u,d for the light-quark current and q — s for 
the strange-quark vector current. A single inter- 
polating field for the nucleon, Sigma and cascade 
octet baryons has the general form 



abcr aT 



[ q r{x)C l5 ql{xM{x), (5) 



where C is the charge conjugation matrix, and qi 
and q 2 are any of the quarks {u,d, s}. For ex- 
ample, to create a proton, we want q\ — u and 
q 2 = d] for the "Er , q± = s and q 2 = d. By cal- 
culating two-point and three-point correlators on 
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Figure 4. Lattice data (circles) for gA, gss and 
g ss and chiral extrapolation (lines and bands). 
The square is the extrapolated value at the phys- 
ical point. 



the lattice with the same baryon operator, we will 
be able to extract the form factors from Eq. [3] 
We solve for Gm,e using singular value decom- 
position (SVD) at each time slice from source to 
sink with data from all momenta with the same 
q 2 and all [i. 

We constrain the fit form to go asymptotically 
to l/Q 4 at large Q 2 and to have G E (Q 2 = 0) = 1: 



Ge 



AQ 2 + 1 



1 



CQ 2 + 1 (1 + Q 2 /Mlf 



(G) 



By trying various combinations of fit constraints 
on our data with squared momentum transfer 
< 2 GeV 2 , we find C is always consistent with 
0; therefore, we set C = 0. The mean-squared 
electric charge radii can be extracted from the 
electric form factor Ge via 



(rl) 



-6)- 



d (G E {Q 2 ) 



Q 2 =Q 



(7) 



'dQ 2 V G E (0) 

In Figure [5] we plot the electric charge radii with 
the neutron and S° omitted, since the vector con- 
served current gives GE,{ n ,s°}(Q 2 ) ~ for these. 
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We see that there is small SU(3) symmetry break- 
ing between the SU(3) partners p and S + (or E~ 
and S - ); their charge radii are consistent within 
statistical errors. Overall, the SU(3) symmetry 
breaking in the charge radii is much smaller than 
what we observed in our study of the axial cou- 
pling constants; for charge radii, the effect is neg- 
ligible. 

We can take a ratio of the electric radii of the 
baryons which coincide in the SU(3) limit; for 
example, p and E + and £~ and S _ . The domi- 
nant meson-loop contribution is suppressed; thus, 
a naive linear fit to a ratio could be a better de- 
scription than fitting individual channels. Follow- 
ing Sec. [31 we use the SU(3) symmetry measure x 
to parametrize the deviation of the ratio from 1 
due to symmetry breaking: 1+X}„=i c nX n , where 
the next-order corrections contribute at the order 
of x N+1 ; taking n = 1, we expect the remain- 
ing effect should be less than 1% in the expan- 
sion. The fit works fairly well for the extrapo- 

By using the ex- 



lation of , ( ffi p and 

VE/E 



r E>S + 



pcrimcntal value of 



and 



we can 



make predictions for (r%)s+ and (r%)s- ■ 0.93(3) 
and 0.501(10) fm 2 respectively. (Using n = 2 
in for the fit yields C2 zero within error and thus 
gives results consistent with the extrapolation us- 
ing n = 1.) 

Studying the momentum-transfer dependence 
of magnetic form factors gives us the magnetic 
moment via 



Ms = G B M {Q 2 = 0) 



(8) 



with natural units 2A ^ , where Mb are the baryon 
(B 6 {N, E, S}) masses. To compare among dif- 
ferent baryons, we convert these natural units 



into nuclear magneton units (In 



2M N 



there- 



fore, we convert the magnetic moments with fac- 
tors of 4t^ 

Mb 

We can obtain the magnetic moments and radii 
from polynomial fitting to the ratio of magnetic 
and electric form factors, Gm /Ge- From the def- 
inition of the electric and magnetic radii, we ex- 
pect that Gm/Ge ~ A + CQ 2 , where the mag- 
netic moment is fx = AGe(0), and C is propor- 
tional to {r\j) — (r%). In the case of n and 5°, 
we use Ge, p and G\ee- m the ratio instead of 
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Figure 5. The electric mean-squared radii in units 
of fm 2 as functions of M 2 (in GeV 2 ) from each 
quark contribution 



Ge,u and G\e,e°- Figure [S] shows the magnetic 
moments of each baryon compared with its SU(3) 
partner: {n,~°}, {£-,»-}. We find 

that as seen in experiment, the SU(3) breakings 
of the magnetic moments are rather small. As 
we go to larger pion masses (that is, as the light 
mass goes to the strange mass), the discrepancy 
gradually goes to zero as SU(3) is restored. But 
even at our lightest pion mass, around 350 MeV, 
the effects of SU(3) symmetry breaking effect can 
be ignored. The fitted results are consistent with 
what we obtained from the dipole extrapolations. 
We examine the radii differences from the quark 
contributions and observe less than 10% discrep- 
ancy. The ratio approach also benefits from can- 
cellation of noise due to the gauge fields, and thus 
it has smaller statistical error. Therefore, we will 
concentrate on the results from this approach for 
the rest of this work. 

We extrapolate the baryon magnetic moments 
using the SU(3) ratios with the SU(3) symme- 
try breaking measure x. Again, the ratio has 
cleaner statistical signal due to the cancellation of 
fluctuations within gauge configurations, and the 
linear extrapolation is not a bad approximation, 
since potential log terms are suppressed. With 
the help of experimental /is- and (i^o |10j. 

we obtain [i p = 2.56(7), ji n = —1.55(8) and 
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Figure 6. Baryon magnetic moments in units of 
Hn as functions of M% (in GeV 2 ). The leftmost 
points are the experimental numbers. 



/is- = —1.00(3). (The fit using up to x 2 terms 
results in a zero-consistent fit parameter ci and 
yields numbers consistent with the above.) 

SU(6) symmetry predicts the ratio /j^p / /i uP 
should be around —1/2. Compared with what 
we obtain in this work, the ratio agrees within 
2(7 for all the pion-mass points. The heaviest two 
pion points have roughly the same magnitude as 
in the quenched calculation 12 . However, at the 
lightest two pion masses, they are consistent with 
the —1/2 value. The difference could be due to 
sea-quark effects, which become larger as the pion 
mass becomes smaller. A naive linear extrapola- 
tion through all the points gives —0.50(10). SU(6) 
symmetry is preserved in the lattice calculations. 

We also check the sum of the magnetic mo- 
ments of the proton and neutron, fi p + fi n , which 
should be about 1 from isospin symmetry. Again, 
the values from different pion masses are consis- 
tent with each other within 2 standard deviations 
and differ from 1 by about the same amount. 
A naive linear extrapolation suggests the sum is 
0.78(13), which is consistent with experiment but 



about 2er away from 1. This symmetry is softly 
broken, possibly due to finite lattice-spacing ef- 
fects. Finer lattice-spacing calculations would be 
needed to confirm this. 

5. Semileptonic Decays 

The hyperons differ from the nucleons by their 
strangeness, and although hyperon decays via the 
weak interaction have been known for more than 
half a century, interest in their study has not de- 
cayed with time. They provide an ideal systems 
in which to study SU(3) flavor symmetry break- 
ing and offer unique opportunities to understand 
baryon structure and decay mechanisms. The 
low-energy contribution to the transition matrix 
elements for hyperon beta decay, B\ — > Bie~V 
can be written in general form as 

M = ^u B2 {Oa+0^)u Bl u el a (l+ l5 )v„. (9) 

From Lorentz symmetry, we expect the matrix 
element composed of any two spin- 1/2 nucleon 
states, B\ and B2, to have the form 



fi(q 2 )l a 



M 2 



M Bl 



o- 



92{q 2 ) p 
<? a p<f 



la 75 



M Bl M Bl 

with transfer momentum q = p B . 2 — p Bl and V, A 
indicating the vector and axial currents respec- 
tively, fi and g\ are the vector and axial form 
factors, and ji and (72 are the weak magnetic and 
induced-pseudoscalar form factors. They are non- 
zero even with B\ — £?2- SU(3) flavor symmetry 
breaking accounts for the non-vanishing induced 
scalar and weak electric form factors, fs and (72 
respectively. Due to the difficulty in disentangling 
experimental form factor contributions, they tend 
to be set to zero. We will be able to determine 
these form factors using theoretical technique and 
will find them to be non-negligible. 

So far, there are only two "quenched" (where 
the fermion masses in the sea sector are infinitely 
heavy) lattice calculations of hyperon beta de- 
cay, and they are in different channels, S — * n 
and S° — > E + . Guadagnoli et al.[T3] calculated 
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the matrix element X — > n with all of the pion 
masses larger than 700 MeV. Sasaki et al. [T3] used 
lighter pion masses in the range 530-650 MeV 
and DWF to look at the S° decay channel. They 
extrapolate the vector form factor fx using the 
parameter 5 — {Mb 2 — Mb 1 )/Mb 2 - In this work, 
we calculate both decay channels and remove the 
quenched approximation, which often causes no- 
toriously large systematic error. Our fermion sea 
sector contains degenerate up and down quarks 
plus the strange. We use pion masses as light as 
350 MeV, which ameliorates some of the uncer- 
tainty in the extrapolation to the physical pion 
mass. To condense our work for this proceeding, 
we will concentrate on the results from £ — > n. 

To obtain V us , we need to extrapolate f\ to 
zero momentum-transfer (we cannot calculate 
this point directly due to the discrete values of 
momentum accessible in a finite volume) and the 
physical pion and kaon masses. Fortunately, f\ is 
protected by the Ademollo-Gatto (AG) theorem 
such that there is no first-order SU(3) breaking. 
Therefore, the quantity deviates from its SU(3) 
value by the order of the symmetry breaking 
term of 0(H' 2 ), where H' is the SU(3) symmetry 
breaking Hamiltonian; the natural candidate an 
observable to track this breaking is the mass split- 
ting between the kaon and pion. Combining with 
momentum extrapolation (using a dipole form in 
this case), we use a single simultaneous fit: 




Figure [7] shows the result from simultaneously 
fitting over all q 2 and mass combinations for 
the X~ — > n decay. The z-direction indicates 
fx, while the x- and y-axes indicate mass and 
transfer momentum. The surface is the fit us- 
ing Eq. [10] with color to indicate different masses. 
The columns are the data and the momentum 
points from different pion masses line up in bands. 
Our preliminary result for fx is —0.95(3) (which 
is consistent with the quenched result [13]: fx = 
-0.988(29)*.*.) 
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Figure 7. Simultaneous extrapolation in q 2 and 
mass 



The axial form factor g\ is not protected by the 
AG theorem. Experimentally, one is interested in 
its ratio with the vector form factor, gx/ fx, at zero 
momentum transfer. Here we adopt naive linear 
mass combinations (M% + M 2 ) and {M\ - M 2 ) 
and momentum dependence to extrapolate and 
find gx(0)/fx(0) = -0.336(52), which is consis- 
tent with the experimental value of —0.340(17). 
Similar extrapolations are applied to other form 
factor ratios, such as / 2 (0)//i(0) = -1.28(19), 
which is consistent with Cabibbo model value of 
-1.297[U]. 

Finally, we turn our discussion toward the 
SU(3)-vanishing weak-electric gi and induced- 
scalar /3 form factors. Figure [8] presents the mo- 
mentum dependence of the ratios 32 (? 2 ) / fx{q 2 ) 
and fz{.q 2 )/fx{q 2 ) for sea-pion masses rang- 
ing 350-700 MeV. The band indicates our 
mass extrapolation in terms of a naive lin- 
ear dependence of M\ — M 2 and M\ + 
M 2 . We find / a (0)//i(0) = -0.17(11) 
and ff 2 (0)//i(0) = -0.29(20), which are 1.5 
standard deviations from zero. Experimen- 
tally, only a combination of axial form factors 
l5i(0)//i(0)-0.133 32 (0)// 1 (0)| is determined. 
They find 0.327(7) (19), which is consistent with 
our result, 0.297(60). 



8 



H.-W. Lin 



0.2 - 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

6 2 (OeV 2 ) 

00 r _ 




-0.8 - 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

(2 ! (GeV 2 ) 

Figure 8. The momentum dependence of the form 
factor ratios /3//1 (top) and 52 //1 (bottom) with 
various pion masses labeled by triangles, squares, 
downward triangles and diamonds from light to 
heavy; the band indicates the extrapolation at 
the physical mass 



6. Summary and Outlook 

The lattice formulation is a powerful method 
for calculating nonperturbative quantities used in 
fundamental tests of QCD. The results presented 
here demonstrate clear progress toward the long- 
term goals of determining the hyperon spectrum, 
structure and decays from first principles. Ongo- 
ing efforts from lattice community using improved 
operators and algorithms (to improve the noise- 
to-signal ratios), finer lattice spacings (to reduce 
systematic error due to discretization), lighter 
pion masses (to reduce uncertainty introduced by 
chiral extrapolation) on dynamical sea quarks will 
continue to better approximate the physical world 
and provide observables that elude experiments. 
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